Drift without flux: Brownian walker with a space dependent diffusion coefficient 



O 

o 
o 

(N 
5h . 

< 



o 



CO 

c3 



I 

C 

O 

o 



> 

o 
o 
in 
o 
o 
o 



-a 
c 

o 
o 



X 
S3 



P. Langon 1 , G. Batrouni 2 , L. Lobry 1 and N. Ostrowsky 1 
1 Laboratoire de Physique de La Matiere Condensee, CNRS UMR 6622, Universite de Nice-Sophia Antipolis, Pare Valrose, 

06108 Nice, France 

2 Institut Non-Lineaire de Nice, CNRS UMR 6618, Universite de Nice-Sophia Antipolis, 1361 route des Lucioles, 06560 

Valbonne, France 

Space dependent diffusion of micrometer sized particles has been directly observed using digital 
video microscopy. The particles were trapped between two nearly parallel walls making their con- 
finement position dependent. Consequently, not only did we measure a diffusion coefficient which 
depended on the particles' position, but also report and explain a new effect: a drift of the parti- 
cles' individual positions in the direction of the diffusion coefficient gradient, in the absence of any 
external force or concentration gradient. 
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Brownian motion of spherical colloidal particles in the 
vicinity of a wall has been extensively studied, both the- 
oretically and experimentally . It has been shown that 
the diffusion coefficients parallel or perpendicular to the 
wall were greatly reduced when the particles were close 
enough to the obstacle, i.e. within distances compara- 
ble to or less than their radius [Q. When the Brown- 
ian particles are trapped in a more confined geometry, 
such as a porous medium, the theory is far more compli- 
cated and few experimental studies have been reported 
in model geometries, where the particles are trapped be- 
tween two parallel walls ^||. In this article, we report 
some new experimental results concerning the Brownian 
motion of particles trapped between two nearly parallel 
walls, so that the confinement, and thus the diffusion co- 
efficient, become space dependent. As a result, we not 
only measure a diffusion coefficient which varies with the 
confinement, but also a drift of the particules' individ- 
ual positions in the direction of the diffusion coefficient 
gradient, in the absence of any external force or con- 
centration gradient. This drift was not accompanied by 
any net particle flux, i.e. statistically the same number 
of particles crossed any imaginary surface in both direc- 
tions. We first discuss the general problem of a Brownian 
walker with a spatially dependent diffusion coefficient to 
explain the origin of the expected drift, and then present 
the experimental set-up and results. 

As in our experiment the diffusion coefficient varies in 
only one direction, say x, we briefly sketch a heuristic 
derivation of the ID Brownian walker algorithm. The 
velocity of a ID Brownian particle subjected to a random 
force and a viscous drag follows the Langevin equation, 



dv(t) 
dt 



= -jv(t)+T(t), 



(1) 



where 7 1 is the velocity relaxation time and T(t) the 
random force per unit mass defined by its mean value 
(r(i)) = and correlation function (T(t)T(t')) = qd(t - 



t'). Using the equipartition theorem it can be shown that 
q is related to the temperature T and the particle's mass, 
m, by the standard relation q = 2-fkT/m. Discretizing 
the random function T(t) over time intervals At >> 7 -1 
allows us to drop in Eq.(l) the inertial term, dv/dt, and 
to replace the velocity v by Ax/ At. Choosing for T(t) 
the simplest random function, T(t) = ±y/q/At, leads to 
the well known Brownian walker algorithm, 



x(t + At) = x(t) ± V2DAi 



(2) 



with D = kT/mj. When the diffusion coefficient D, 
i.e. when the temperature T and/or the drag coeffi- 
cient 7 become position dependent, the above algorithm 
needs to be clarified. During each time interval At, the 
walker makes a step to the right or to the left, but should 
the length of this position dependent step, \J2D At, be 
computed at the departure point x(t) = x, the arrival 
point x(t + At) = x + Ax or at any point in between? 
These mathematical choices, often referred to as the 
Ito/Stratonovitch conventions Q, model different phys- 
ical situations and the choice of convention is dictated 
only by the physics. We denote by D(x + a Ax) the dif- 
fusion coefficient appearing in Eq.(2) where a = 0,1/2 
and 1 correspond to the Ito, Stratonovitch and isother- 
mal choices respectively. As we will show, this last case 
models a situation where the temperature, T, is uni- 
form but the drag coefficient, 7, is space dependent. 
Using in Eq.(2) the limited expansion, D{x + aAx) 
D(x) + a(dD/dx)Ax with Ax = ±^2D(x)At, yields 
the algorithm for a Brownian walker with a position de- 
pendent diffusion coefficient: 



x(t + At) 



dD 



x(t) ± J2D(x(t))At + a — At. 

dx 



(3) 



Depending on the value of a, this model has very different 
implications concerning the equilibrium distribution of 
the Brownian walkers, their individual drift x(0))) 
and their net flux. 
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Averaging Eq. (3) over a large number of walkers shows 
the average position of a Brownian walker is no longer 
zero since the diffusion gradient term acts as an external 
force leading to particle drift. If this gradient is assumed 
to be constant, this drift increases linearly with time as 



(x(t) - x(0)) = a^-At. 

dx 



(4) 



A first intuitive, but misleading, idea would be to con- 
clude that the particles will migrate in the direction of 
the diffusion gradient, leading, therefore, to a concen- 
tration gradient. This is actually incorrect as we now 
show. Starting with a uniform particle distribution po, 
we check if this corresponds to an equilibrium state by 
determining the particle flux through an imaginary sur- 
face, S, placed perpendicular to the diffusion coefficient 
gradient at coordinate x (see Fig. 1). During a time in- 
terval At, all particles crossing S from the left (or right) 
are half of those included in the volume SL r i g ht (SLi e f t ) 
where L r i g } lt (Li e f t ) is the right (left) step terminating 
at x taken by a walker during At. The net particle flux 
to the right will thus be 



J 



Po SL 



right 



SL 



left 



SAt 



(•5) 



Equation (2) allows computing the length of these two 
steps which both end at the same point x: 



dD 



LrigM = J2D(x)At± (a - I)— At, 

iijt ' ax 



leading to the particle flux 



J 



-p (l - a) 



dD_ 

dx 



(0) 



(7) 
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FIG. 1. Particle flux for Brownian walkers with a step 
length depending on arrival position (a — 1) or departure 
position (a = 0). 

As a result, in the situation of maximum drift where 
a = 1, this flux will vanish (see left part of Fig. 1), mean- 
ing that the uniform particle distribution corresponds to 
an equilibrium. According to Boltzmann, this should 



correspond to an isothermal situation, the diffusion co- 
efficient gradient arising only from a pure hydrodynamic 
effect, the spatial dependence of the drag coefficient 7. 
For all the other values of a, the flux will be negative, 
leading to a concentration gradient of the particles in the 
direction opposite to that of the diffusion coefficient gra- 
dient. The maximum flux (and zero drift) is obtained 
for a — 0, as shown on the right part of figure 1. Equa- 
tion (7) may be generalized to the case where the particle 
distribution is position dependent: 



dx dx 



(8) 



which leads to the well known generalization of Fick's 
law for a space dependent diffusion coefficient 

To confirm these results, we performed simulations of 
Brownian walkers following algorithm (3). We found that 
only the a = 1 case leads to a uniform distribution of par- 
ticles with no net flux through any given surface while, 
at the same time, the average individual positions exhibit 
a drift in the direction of the diffusion coefficient gradi- 
ent according to Eq(4). This situation of "drift without 
flux" may be compared to the equilibrium situation of 
Brownian particles subjected to an external force, such as 
their weight: If one follows the motion of individual par- 
ticles, an average downwards drift is observed; however, 
there is no net flux because of the vertical concentration 
gradient. In our isothermal case, the drift of individual 
particles from lower to larger D(x) region does not lead 
to a net flux because particles in the larger D(x) region 
diffuse further than particles in the lower D(x) region. 
This physical situation imposes the choice of a = 1 in 
algorithm (3), so that a particle coming from a low D(x) 
region makes a right step just equal to the left step of that 
particle coming from a high D(x) region and arriving at 
the same point (see left part of Fig. I). 
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FIG. 2. Experimental set-up. Rectangular inset is an 
enlarged top view showing the center of the cell with the cir- 
cular excluded volume and the observation frame. Round 
inset explains the two contributions to the change in diffusion 
coefficient when a particle moves a distance dx. 
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We have set up an experiment to check this drift with- 
out flux prediction where particles, observed under a mi- 
croscope, undergo Brownian motion in a confined geom- 
etry. When the confinement, e, is of the order of the 
particle size, a, the diffusion coefficient strongly depends 
on the value of e/a. As the confinement was position de- 
pendent, we observed Brownian walkers in an isothermal 
situation but with a spatially dependent D(x) due to a 
purely hydrodynamic effect. 

Polystyrene spheres, of radius a = 1/im, were sus- 
pended in a mixture of H2O + D2O so as to cancel any 
sedimentation effects. Addition of a surfactant (2.2g/l 
of SDS) helped minimize particle aggregation or adhe- 
sion to the walls. A drop of this mixture was placed 
between a flat disk and a planar convex lens, Fig. (2), 
of curvature radius R = 15.5mm separated by an clastic 
O-ring. The flat and convex surfaces were then brought 
into contact at the center of the cell by gently squeezing 
the elastic joint, the remaining air providing the neces- 
sary sample compressibility. The spacing, e, between the 
flat and curved wall depends on the distance r from the 
center of the cell as e = r 2 /2R. The contact between 
the two walls as well as the dependence of e on r were 
carefully measured by monitoring the Newton rings ob- 
served under the microscope. We used as a light source 
a new super-radiant diode B whose coherence length is 
less than 100/im. This was important as this coherence 
length was long enough to observe the desired Newton 
rings, but short enough to avoid any other interference 
patterns due to all the cell interfaces, which were visi- 
ble with an ordinary diode laser and which completely 
masked the relevant signal. The horizontal Brownian 
motion of the polystyrene balls was observed through a 
microscope equipped with a long range objective of mag- 
nification 50A, followed by a CCD camera coupled to 
the microscope via an eyepiece of magnification 8X. The 
video signal was processed in real time by a computer, 
which recorded, every 3s the horizontal position, size and 
shape of all objects in a rectangular frame 65/imA 100/im. 
This time interval was long enough to allow the image 
analysis of all particles present in the frame and small 
enough so that the particle's average displacement was 
only a fraction of their diameter. 

As the particles' confinement, e, was related to their 
distance, r, from the center of the cell, we were able to ex- 
plore different confinement regions by moving the obser- 
vation frame in the horizontal plane. The explored e var- 
ied from 2.5/im to 11/im. The vertical Brownian motion 
of the particles over this small vertical range could not be 
monitored. However, we took that motion into account 
when interpreting the data by averaging the particle's 
vertical position over the confinement range. The vol- 
ume fraction of polystyrene balls, of the order of 1%, was 
chosen so as to allow the monitoring of a fairly large num- 
ber of particles at the same time (around 30 for e = 3/im) 
to improve the statistics in the data analysis. Following 



the particles' positions from one frame to another, the 
program analyzed a great number of trajectories (more 
than 10 5 for each run). When two particles got closer 
than twice their diameter, the program treated them as 
"dimers", their trajectories as "monomers" were ended 
at that time and were no longer used to determine the 
diffusion coefficient or the drift. In that way, the role of 
particle interactions, which have a range smaller than a 
couple of particle diameters, could be safely ignored. 
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FIG. 3. Du/Dq with respect to relative confinement e/2a 
. Open squares are the experimental data, black dots were 
calculated by the collocation method, and dotted and solid 
lines follow analytical approximate solutions. 

The confinement dependence of the diffusion coefficient 
was determined as follows. A given observation frame 
was divided into 3 zones (see inset, Fig. (2)), each corre- 
sponding to an approximately constant e. For each zone 
we averaged the particle's displacement squared, either 
in the x, or in the y directions, as a function of time, and 
checked that indeed they followed the usual diffusion law, 



= (y 2 ) =2D||(e)t 



(9) 



By moving the frame to different locations, we were 
able to explore a range extending from e/2a = 1.2 to 
11. The bulk diffusion constant, Do, was determined 
using the well-known relation Dq = kT/Gnria, where 
1] = 0.99A10 _3 SI is the viscosity of the water- heavy 
water mixture, yielding Dq = 1.92X 10~ 13 m 2 /s. The ex- 
perimental values for D\\ / Dq are shown in Fig. (3) (white 
squares) and fit remarkably well the available theoretical 
predictions (black dots) using the collocation method (?j 
averaged over all the possible vertical positions z of the 
particle for a given e, i.e. with a < z < (e — a). For 
comparison, we also plotted (solid line) the analytical 
solution obtained using the Faxen expression || for the 
position dependent drag of a particle moving parallel to 
a single wall, then adding the effect of each of the two 
walls, and averaging over the vertical position z. This 
solution clearly overestimates the reduction of the diffu- 
sion coefficient of a particle trapped between two parallel 
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walls, particularly as the relative confinement, e = e/2a 
reaches its lower limit 1. We also plotted (dashed lines) 
in Fig. (3) the analytical expression for the drag of a par- 
ticle trapped just in the middle of two parallel planes. ||] 
The impossibility of averaging over z an expression only 
known for z — e/2 explains the observed discrepancy 
which goes to zero as the relative confinement approaches 
its limit, 1, where the z-average becomes irrelevant. 

To demonstrate the existence of an individual drift of 
the particles, we fixed the center of the observation frame 
at a position y = and x — 300^m, corresponding to an 
average e/2a — 1.5 so that all particles present in the 
frame were outside the excluded volume (i.e. e > 2a), 
and had a diffusion coefficient with the largest x depen- 
dence, but no y dependence (to first order). For the 
determination of (x(t) — x(0)) and (y(t) — 2/(0)), each 
trajectory was divided into independent paths lasting a 
time t, each contributing to the evaluation of the aver- 
age drift during time t. The results are shown in Fig. 
(4), and reveal a drift in the Brownian walker position 
along the x direction, and none in the y direction along 
which the diffusion coefficient may be considered as con- 
stant. The statistics of the results clearly deteriorates as 
time increases: After recording trajectories for typically 
a dozen hours, more than a hundred thousand indepen- 
dent segments contributed to the determination of the 
drift at short times, whereas only up to a few thousand 
independent segments were left for t = 200s. This is due 
to the fairly high particle concentration which lowers the 
lifetime of a "monomer" (time during which a particle 
doesn't approach another one to within 2 particle diam- 
eters), but which was chosen as a compromise to have 
good statistics at short times while allowing us to follow 
each particle during a reasonable time, fixed at 200s. It 
should be pointed out that in order to avoid any bias in 
the statistics, for a trajectory segment to be valid and 
included in the statistics, the position of a walker at in- 
stant t = had to be inside a region 15/im away from the 
edges of the observation frame. This condition ensured 
that after diffusing for 200s, the walker had less than 
0.5% chance to have covered 15/xm, and was thus still 
present in the observation frame. Failure to impose this 
condition resulted in the observation of a spurious drift, 
in the opposite direction, due to an artificial selection of 
walkers because of the experimental boundary conditions 
(limits of the observation frame) . 

To compare our experimental results with the theo- 
retical predictions, Eq. (4), we evaluated the diffusion 
gradient encountered by the walkers present in the obser- 
vation frame. It is important to realize that as a walker 
moves a distance dx (see round inset in Fig. 2), its dif- 
fusion coefficient Du varies first because the confinement 
varies (path (1) parallel to the bottom wall, i.e. at con- 
stant z), and second because at constant confinement, 
the particle's height z changes (path (2)). Adding both 
contributions and averaging over the vertical position z 



of the walker yields: 

, dD\\ (e, z) . x ,dD\\(e,z) x ,dD\\(e,z) 
' dx ' z _ d~e ' z + 2R^ dz 

(10) 

Using our experimental data and results of the colloca- 
tion method, we found (dDn(e, z)/dx) z « 2.2A10 _9 m/s. 
This value of the slope is used to plot the straight dotted 
line on Fig. 4. The experimental data are thus in good 
agreement with the predicted drift corresponding to the 
expected a = 1 value. 
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FIG. 4. Average position of the walkers as a function of 
time, along (black dots) and perpendicular (open squares) to 
the diffusion gradient. 

Finally, to claim drift without flux, it is not sufficient 
only to demonstrate the drift: we must also demonstrate 
the absence of flux. If a flux was due to the observed 
drift dD/dx, we would expect a radially outwards flux 
of pdD/dx particles, which would empty our observation 
screen in less than a day. Furthermore, if this flux were to 
be balanced by a concentration gradient, one can show 
that a concentration change by 30% over a distance of 
60/im would be necessary. Experimentally, we observed 
no flux and no concentration gradient over a period of 
a week or more, which is consistent with the Boltzmann 
requirement of a uniform concentration in the absence of 
a temperature gradient. 
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